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Abstract 

The relation between dark matter halos and the loci of star formation at high redshift is a pressing question in contemporary cosmology. 
Matching the abundance of halos to the abundance of infrared (IR) galaxies, we explore the link between dark matter halo mass (M h ), 
stellar mass (M*) and star-formation rate (SFR) up to a redshift of 2. Our findings are five-fold. First, we find a strong evolution of 
the relation between M t and SFR as a function of redshift with an increase of sS FR = S FR/M* by a factor ~30 between z=0 and 
z= 2.3. Second, we observe a decrease of sSFR with stellar mass. These results reproduce observed trends at redshift z>0.3. Third, 
we find that the star formation is most efficient in dark matter halos with M h = 5x10" M G , with hints of an increase of this mass with 
redshift. Fourth, we find that SFR/M/, increases by a factor ~15 between z, = and z = 2.3. Finally we find that the SFR density 
is dominated by halo masses close to ~7xl0" M Q at all redshift, with a rapid decrease at lower and higher halo masses. Despite its 
simplicity, our novel use of IR observations unveils some characteristic mass-scales governing star formation at high redshift. 
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1. Introduction 

The loci of star formation at high redshift is one of the salient 
questions of galaxy formation theory. In particular, uncertainties 
in the relation between the host halo mass and the IR luminosity 
function (LFIR) currently limit models of the clustering of 
star-forming galaxies as revealed e.g. by the anisotropics of the 
cosmic infrared background (CIB). This comes from the fact 
that the clustering of galaxies is directly driven by the mass of 
their host halos. The so-called abundance matching technique 
is another valuable tool that provides insights on this relation 
( |Vale & O striker 2004). While unexplored in the IR, a recent 
improvement of LFIR measurements makes this approach now 
both promising and timely. We will illustrate its power in this 
letter. In particular, working in the IR provides the unique 
possibility to constrain the star-formation rate (SFR) directly 
without any extra assumption regarding the physical processes 
that drive the star formation. 

Considering the latest observational and theoretical develop- 
ments in the IR, we introduce in Sect. [2] the abundance match- 
ing technique. We describe in detail the mass and luminosity 
functions used in Sect. [3] Matching abundances allows us to in- 
fer relations between halo and stellar mass (Sect. |4.1| i, SFR and 
stellar mass (Sect. |4.2| i and SFR and halo mass (Sect. |4.3| l. In 
these sections, we will illustrate in particular how this simple but 
well-informed approach satisfyingly reproduces trends observed 
in the IR, but also at other wavelengths. Focusing on what is re- 
quired for the modeling of the CIB, we exclude the so-called 
starbursting galaxies here because they are non-significant out- 
liers. Critical in the context of abundance matching, this assump- 
tion is well motivated and discussed in Sect.[3j 
Throughout this paper, we adopt a Chabrier (2003) initial mass 



function (I MFfj and a ACDM WMAP-7 cosmology ( |Larson 
|etal.|2011| ). 

2. Simple abundance matching 

Our goal is to find an observationally supported relation between 
the SFR of galaxies, their stellar and their host halo mass. While 
this relation is statistical in nature, the abundance matching tech- 
nique (e.g. Conroy & Wechsler 2009) provides a way to relate 
the means of this relation. Observations suggest a small scatter 
around this relation and we will therefore neglect it: 0.16 dex for 
the relationship between stellar and halo mass (More et al. 2009 ), 
and 0.2 dex between the stellar mass and the SFR (Rodig hiero] 
et al. (201 1 1; see also Behroozi et al. ( 2010| i for a more general 
discussion of scatter). 

The key assumption of this matching technique is the exis- 
tence of a continuous and monotonic relation between the halo 
mass, the stellar mass, and the SFR: SFR = /(M*), SFR = 
g(M h ) and M* = h{M h ), where SFR is the SFR, M* the galaxy 
stellar mass (SM), and Mh the dark-matter halo mass (HM). 

We solve for the functions /, g and h that satisfy the follow- 
ing relations between the infrared luminosity function, the halo 
mass function, and the stellar mass function: 

n LlR (> L m = K 1 x /(M*)) = ««,(> M*) (1) 
n Lm (> L m =K~ l xg (M/,)) = n Mh {> M h ) (2) 
%,(> M* = h(M h )) = n Mh (> M h ), (3) 

where ni IR is the number of galaxies per unit comoving vol- 
ume brighter than the bolometric infrared luminosity Lm, and 



Stellar masses and SFR computed assuming a Salpeter ( 1955 1 IMF 
can be converted to Chabrier (2003) IMF applying a -0.24 dex correc- 
tion. 
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and riM h are the number of galaxies more massive than 
the stellar mass and the halo mass M/„ respectiveljj^] K is 
the |Kennicutt| ( |1998| l factor, which links the infrared bolomet- 
ric luminosity integrated between 8 and 1000 /mi and the SFR: 
As we discuss below, a key assumption 
in this matching is that all (sub)halos contribute to the IR. 



3. Mass and luminosity functions 

The key inputs to our discussions are the bolometric infrared 
luminosity function (LFIR), the stellar mass function (SMF) and 
the halo mass function (HMF). 

It is well known that the LFIR is inconsistent with the 
Schechter (|1976|l form an d presents an excess at the bright end 
(e.g. Saunders et al.| 1990 1. However, it has been shown that this 
excess is mainly caused by starbursting galaxies, which are well 
outside the SFR-SM main sequence. |Yun et aT ( 2001| l showed 
that the 60 /mi LF of the IRAS 2 Jy sample can be fitted by the 
sum of two Schechter functions. More recently, Sargent et al. 
(2011 ) showed that the LFIR can be decomposed as a sum of 
two approximate Schechter functions. The first one is respon- 
sible for ~85% of the IR energy output, and is caused by the 
so-called main-sequence objects that have very similar specific 
sSFR at a given M* and redshift. The second one dominates the 
bright-end of the LF, and is caused by starbursting galaxies with 
an sSSFR about 0.5 dex higher than the main sequence. These 
starbursts contribute only ~2% to the mass-selected star-forming 
galaxies number density ( Rodighiero et al.|201 1} , and ~15% to 
the SFR density (Sargent et al. 2011), and consequently to the 
CIB. The former guarantees that starbursts constitute a negli- 
gible fraction of IR objects, i.e. do not contribute significantly 
to the IR abundances. The latter suggests that they can be ne- 
glected for our purpose. Additionally, because a bimodal rela- 
tion between LFIR and sS FR would violate our hypothesis of a 
monotonic relation between masses and SFRs, we henceforth ig- 
nore the starburst components. Focusing on the main-sequence 
population, we therefore write the LFIR as a Schechter ( |1976| l 
function 



dN 
dL 1K 



^1R 



exp 



1 



(4) 



where dN/dLm is the LFIR, and 0*(z) and L*(z) are the charac- 
teristic density and luminosity at a given redshift, respectively. 
We use the evolution of 0*(z) and L*(z) measured by Magnelli 
etaLlpOTT) up to z=2.3. 

In contrast to the LFIR, the SMF is much closer to a 
|Schechter| ([1976} function. |Perez-Gonzalez et al.] ( |2008| l have 
measured the SMF up to z=4 and we used their parameters of 
the Schechter function fit. We performed our abundance match- 
ing at the same redshift as the center of the redshift bins of|Perez-| 



Gonzalez et al.|(|2008]l. We checked that there is no inconsistency 



between the SFR derived from our LFIR and that derived from 
the SMF following [Wilkins et al |(j2Q08>. 

For the HMF, we used the |Tinker et al. ( 2008 1 z-dependent 



HMF and the associated sub-halo mass function proposed in 
|Tinker & Wetzer] ( pOlO) with A = 200. Our results are fairly 
insensitive to the inclusion of the sub-halo mass function. 



2 Note that n can be deduced from the luminosity/mass function: 

n x {> X mi „) = l^Jv dX > where X is L '«' M *' or M h- 

3 The Kennicutt factor was re-scaled from the Salpeter to the 
Chabrier IMF. 



4. Results 

4.1. Relation between stellar mass and halo mass 

The abundance matching between the HMF and SMF has been 
extensively studied for z<l (e.g. |Conroy & Wechsler|2009|l and 
has been discussed in detail up to z~4 by |Behroozi et al.| l |2010) . 
Because M*/M/, as a function of M/, is one basic element of 
our analysis, we first investigated it using our mass functions. 



the measurements at z~0 and z~l ( 


Mandelbaum et al.|2006; 


|Conroy et al.|2007t|Reyes et al.|201 1 


I. However, constraints are 



poor except at z = for M h > 10 u M . We find that M+/M/, 
peaks near M/,=10 12 M Q and decreases at lower- and higher- 
mass. This suggests that the integrated star-formation efficiency 
(= M+/Mf,) decreases for smaller as well as higher-mass halos. 
Fig. [T] also indicates that the mass corresponding to the peak 
of integrated star-formation efficiency evolves with redshift, be- 
ing a factor of 8 larger at z=2.25 than at z=0. We extended our 
analysis to higher redshifts using the latest SMF measurements 
( |Caputi et al.|201 l[|Kajisawa et al. 200 9]) and con firm this trend. 



However, as discussed in |Behroozi et ak] ( |2010| l, these results 
have to be taken with care because the current uncertainties make 
these evolutions quite uncertain. 



O 




Figure 1. Ratio between stellar and halo mass as a function of 
halo mass from z=0 to z=2.25. Solid lines: our results based on 
abundance matching technique. Triangles: measurements from 
galaxy-galaxy lensing (Mandelba um et al.|2006|l. S quares: mea- 
surements from satellite galaxies ( Conroy et al. 2007 1 at z~0 
(black) and z~l (green). Diamonds: measurements from com- 
bined kinematics and galaxy-galaxy lensing (Reyes et al. 201 l| l. 
Typical uncertainties on the results of abundance matching are 
0.2 dex. 



4.2. Relation between star-formation rate and stellar mass 

Matching the abundances derived from the LFIR and the SMF, 
we plot in Fig.|2]the sSFR as a function of M* and redshift. We 
observe a strong increase of the sSFR with redshift (about 1 .5 or- 
der of magnitude between z = and z = 2), which qualitatively 
agrees with observations (El baz et al.|201 1| Karim et al.|2011 



We also see a slight decrease of the sSFR with M* (0.5 dex 
between M* = 10 9 ' 5 M Q and M* = 1O 1L5 M ), which again 
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agree s with the observations (Rodighier o et al.|201 l[|Karim et al.| 



201 1 1. To further illustrate this qualitative agreement, we also 
plot in Fig. [2] the most recent sSFR measurements. We repro- 
duce their trends quite well. We note the large dispersion be- 
tween the measurements, however. For example, |Karim et al. 
( 201 1 J and Rodighiero et al. ( 201 1] > differ significantly at high- 
mass at z ~0.3, and at all masses for z~l (note however that 
the redshift bins are quite different). Also, at higher redshift the 
measurements seem to differ by factors ~2 around 10 M . 
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Figure 2. sSFR as a function of stellar mass at various red- 
shifts. Solid lines: our results based on the abundance match- 
ing technique (we use dashed lines at z < 0.3 where neglecting 
the quiescent galaxies have a significant impact, see Sect. B|. 
Filled squares: measurements from |Karim et al.| ( |201 1) obtained 
by radio stacking at z =0.3, 0.5, 0.7 0.9, 1.1, 1.4, 1.8 and 2 .25. 
Filled diamonds: measurement from Rodighiero et al. ( 2010 1 ob- 
tained by Herschel/PACS stacking at z =0.25, 0.75 , 1.25, and 
1.75. Filled circles: measurements from Bauer et al. (201l| ob- 
tained using UV measurements at z =1.75 and 2.25. Filled tri- 
angles: measurements of Oliver et al. (2010) from Spitzer/MlPS 
stacking at z =0.1, 0.3, 0.5, 0.7, 0.9, 1.1, 1.35, 1.75. The er- 
ror bars have not been plotted for clarity, but they roughly 
match the symbol size. Thin-dashed and dotted lines are the lim- 
its between the normal and LIRG (luminous infrared galaxies, 
10 11 L < L IR < 10 12 Lq), and the LIRG and ULIRG (ultra lu- 
minous infrared galaxies, Ljr > 1O 12 L ) regime, respectively. 



4.3. Relation between star-formation rate and halo mass 



We have seen in Sect. |4. 1| and |4.2| that matching abundances 
leads to elucidating connections between the stellar and halo 
mass on the one hand, and the stellar mass and the sSFR on the 
other hand. It is therefore natural to investigate the connection 
between the halo mass and the SFR. We can make this connec- 
tion in two ways. In Fig. [3] (upper panel), we plot the redshift 
evolution of the ratio between SFR and halo mass as a function 
of halo mass. The solid lines use a direct matching between 
the LFIR and the HMF The dotted-dashed lines stem from a 
combination between the function linking the stellar mass and 
the halo mass, and the power-law fit of the sSFR-M* relation of 
|Karim et aLlpOTT) (SF R(M h ) = f(h{M h )) with / taken from 
the Karim et al. ( 20 1 1| > fit and h from the results presented in 



Sect. |4.1[ ). In these two cases, we find a very strong evolution of 
the ratio between the SFR and the halo mass with redshift (1.2 
order of magnitude between z=0 and z=2.25). We also see an 
increase of the halo mass where the instantaneous star formation 
is the most efficient, from M/^xlO 11 M G at z=0 to lxlO 12 M G 
at z=2.25. This evolution is even more pronounced with the 
second method. The relationship between Mi, and SFR/M/, is 
close to a log-normal distribution with a typical width ~0.75 
(~0.65 if we use the second method), but presents an excess 
at high mass (M/, > 10 13 M Q ). This differs notably from the 
SFR-M/, relation found by |Conroy & Wechsler| ( |2009] > at high 
mass (M/, > 2 x 10 12 M o ). However, at these high masses the 
conversion between their stellar mass growth and SFR is less 
reliable (see their Sect. 3.2). In contrast to |Conroy & Wechsler 
(2009) we made no assumption regarding the relation between 
star formation and halo growth. 




10 10 
M„ [M e ] 

Figure 3. Upper panel: Ratio between SFR and halo mass as 
a function of halo mass at various redshifts. Solid lines: results 
from a direct abundance matching of the LFIR and the HMF 
(we use dashed lines at z < 0.3 where neglecting the quiescent 
galaxies has a significant impact, see Sect. |5j. Dotted-dashed 
lines: results from the abundance matching between the SMF 
and the HMF, combined with the best-fit of the sSFR-M* rela- 
tion by |Karim et al.| ( [201 1) . Thin-dashed and dotted lines are the 
limits between the normal and LIRG, and the LIRG and ULIRG 
regime, respectively. Lower panel: contribution of the various 
halo masses to the SFR density. 



From the SFR-M/, relation, we can easily derive the differ- 
ential contribution of each halo mass to the global SFR density 

PSFR- 



dp, 



'SFR 



d 2 N 



dlog(M h ) dM h dV 



xSFR(M h ), 



(5) 
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where d 2 N/(dM h dV) is the HMF and SFR(M h ) is the SFR as- 
sociated to a given halo mass provided by the abundance match- 
ing. Fig. 3]shows the evolution of this observable with halo mass 
and redshift. We find that the SFR density is dominated by halo 
masses close to ~7xlO n Mo, in agreement with the analysis 



very crude approximation of a suppression of acc retion and SFR 



in halo with M h < 10' 1 M Q , Bouche et al. (2010) succeeded in 



of Conroy & Wechsler 



(20091 (their Fig. 10). The distribution is 



quite sharp with a rapid decrease at smaller and higher mass. We 
see specifically from this figure that the bulk of star formation at 
z<2.3 never occurs in small systems. 

5. Discussion 

We presented a first application of the abundance matching tech- 
nique to current infrared data. Our exercise turned out to be sur- 
prisingly elucidating. Despite its simplicity, this technique pro- 
vides a reasonable picture of the evolution of the SFR with the 
stellar mass, in qualitative agreement with the multi-wavelength 



reproducing the mass and redshift dependences of the SFR-M* 
and Tully-Fisher relations from z~2 to the present. The exact 
physical processes leading to this mass floor are still under de- 
bate. Again, with our simple abundance matching approach, the 
drop of SFR in low-mass halos appears naturally. These charac- 
teristic mass-scales also roughly agree with those required by the 
constraints coming from the clustering modeling of star forming 
galaxie s (e.g. Cooray et al. ( 2010 >; Planck Collaboration XVIII 
( 201 1\ ) although the detailed relation between our characteristic 
masses and theirs require additional investigations. 

Because the halo mass is the prime quantity that defines 
the matter clustering on large scales, the link between SFR and 
halo mass is not only fundamental for understanding the evo- 
lution of the galaxies, but also the CIB anisotropies. The cur- 
rent CIB anisotropy models, based on halo occupation distribu- 



measurements at z>0.3 (Sect.H~2]i. It confirms the strong link be- tion prescriptions (|Penin et al.|201 1 ||Planck Collaboration XVIII 



tween these two quantities. The success of this analysis relies on 
the hypothesis that the contribution of main-sequence outsider 
objects to the IR energy output is negligible. Matching the LFIR 
and the SMF including the starburst contribution into LFIR leads 
to a very strong increase of the sSFR-M* relation for M* > 10 11 
M , which is not observed. Similarly, it artificially increases the 
sSFR-M/, relation for M h > 3 x 10 12 M . The tightness of the 
sSFR-M* relation indicates that the SFR is not driven by merger- 
induced starbursts but instead by a continuous mass-dependent 
process that is gradually declining with time. Reproducing the 
evolution of this relation is a challenge for models of galaxy for- 
mation (e.g. [DavS]2008] |Damen[etar]2009]l. 

Our analysis neglects the presence of quiescent massive 
galaxies at high mass. By comparing the matched abundances 



with or without quiescent galaxies in the SMF using the Ilbert 
|et al.| ( |2010[ ) measurements, we found that their contribution is 
negligible (< 0.2 dex) for z > 0.3, but reach 0.25 dex at z = 0.3. 
A smaller effect (~0.15dex) at low z is found using the mea- 
surements of jvan der Wei et aL] ( |2010[ ). Because the CIB is dom- 
inated by z > 0.3 galaxies, neglecting the contribution of quies- 
cent galaxies has little impact on its modeling. 

Fig. [3] demonstrates that the halo mass with the most effi- 
cient star formation does not evolve strongly with redshift, with 
a decrease of only a factor ~3 since z^2.3. This may contradict 
with one meaning of the downsizing phenomenon, which implies 
which stars are being formed in preferentially smaller systems 
at later times. This lack of evidence has also been observed by 
|Conroy & Wechsler] ( |2009| > up to z~l. This effect could be pro- 
duced by a poor estimate of the faint-end slope of the LFIR. To 
investigate this, we repeated our analysis using higher slopes (up 
to or=1.8) and found that our conclusion still holds. 

The existence of a characteristic mass-scale at high-mass 
{Mi, a 10 12 M Q ) has been known since the late 1970s and early 
1980s (e.g. |Rees & Ostriker|[T977l [SHk|[T977l |White & Rees] 
1978). More recently, cosmological simulations also suggest the 



existence of a critical halo mass, with an important role in shap- 
ing up the main properties of galaxy. Below this mass, galaxies 
are built by cold flows associated with efficient star formation. 
Above this mass, cooling and star formation are shut down by 
shock-heating triggered feedback (e.g. Dekel & Birnboim 2006). 
Our abundance matching naturally finds this high-mass cutoff 
for the star-formation (although it is quite a smooth transition). 
At the low-mass end, there is also strong evidence for a drop 
in the efficiency of galaxy formation (e.g. Shankar et al.| [2006 1 
van den Bosch et al.||2007] |Baldry et al.||2008||Kravtsov||201o[ 



Moster et al.|2010||Guo et al.|2010[|B"ou che et a l.|2010| i. With a 



|2011| Amblar d~et al.||2011| >, do not take into account that the 
brighter infrared sources tend to be in more massive halos (al- 
though Shang et al. |20 1 1 1 attempted to address this). Our analy- 
sis showed that the relation between SFR and the halo mass has 
a fairly constant shape from z=0 to z=2, but significantly shifts 
with the redshift, the characteristic halo mass varying by a fac- 
tor ~3 with redshift (see Fig. [3}. With improved models of CIB 
anisotropies, we can hope to accurately constrain this relation 
and its evolution with halo mass and redshift. This modeling may 
be able to decrease the tension between wavelengths discussed 
in the |Planck Collaboration XVIII| ( |20TT) . 
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